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ABSTRACT 

We investigate inertial mode oscillations of slowly and uniformly rotating, 
isentropic, Newtonian stars. Inertial mode oscillations are induced by the Coriolis force 
due to the star's rotation, and their characteristic frequencies are comparable with the 
rotation frequency $7 of the star. So called r-mode oscillations form a sub-class of the 
inertial modes. In this paper, we use the term "r-modes" to denote the inertial modes 
for which the toroidal motion dominates the spheroidal motion, and the term "inertial 
modes" to denote the inertial modes for which the toroidal and spheroidal motions have 
comparable amplitude to each other. Using the slow rotation approximation consistent 
up to the order of 0^, we study the properties of the inertial modes and r-modes, 
by taking account of the effect of the rotational deformation of the equilibrium on 
the eigenfrequencies and eigenfunctions. The eigenfrequencies of the r-modes and 
inertial modes calculated in this paper are in excellent agreement with those obtained 
by Lindblom et al (1999) and Lockitch &; Friedman (1998). We also estimate the 
dissipation timescales due to the gravitational radiation and several viscous processes 
for polytropic neutron star models. We find that for the inertial modes, the mass 
multipole gravitational radiation dominates the current multipole radiation, which 
is dominating in the case of the r-modes. It is also found that the inertial mode 
instability is more unstable than previously reported by Lockitch & Friedman (1998), 
and survives the viscous damping processes relevant in neutron stars. 

Subject headings: instabilities — stars: neutron — stars: oscillations — stars: rotation 

1. Introduction 

Since Andersson (1998) suggested and Friedman & Morsink (1998) analytically verified that 
the r-modes of rotating stars are unstable due to the gravitational radiation reaction, much 
attention has been paid to oscillation modes of rotating neutron stars because of their possible 
importance in a field of astrophysics. Lindblom et al. (1998) argued that due to this instability, 
the maximum angular rotation velocity of hot young neutron stars is strongly restricted. Owen et 
al. (1998) also suggested that the gravitational radiation emitted from hot young neutron stars 
due to the r-mode instability is expected to be one of the potential sources for the gravitational 
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wave detectors, e.g. LIGO. Following these earlier investigations on the r-mode instability, there 

have appeared a number of papers which study the properties of the r-modes of rotating stars 
(e.g., Kojima 1998, Anderson et al 1998, Kokkotas & Stergioulas 1998, Lindblom k, Isper 1998, 
Beyer & Kokkotas 1999, Lindblom et al 1999, Kojima & Hosonuma 1999). 

It is well known that the r-modes form a sub-class of the modes called "inertial modes" , for 
which the restoring force is the Coriolis force and the characteristic frequency is comparable to the 
angular rotation frequency $7 of the star (e.g., Greenspan 1964). As the prominent characteristics 
of the r-modes, we know that the toroidal motion dominates the spheroidal motion in the 
velocity field, and that the frequency u> observed in the corotating frame of the star is given by 
u = 2m^/{l{l + 1)) in the limit of f2 — > 0, where m and I are the indices of a spherical harmonic 
function Yj^ representing the toroidal component of the velocity field (e.g., Papaloizou & Pringle 
1978 , Provost et al 1981, Saio 1982). As shown by Lee et al (1992), however, the inertial modes 
in general have the spheroidal component of the velocity field comparable with the toroidal 
component and do not have any analytical formula to give w/O in the limit of O — ^ 0, except for 
the case of the Maclaurin spheroids (Lindblom & Ipser 1998; see also Bryan 1889). In this paper, 
to distinguish these two classes of inertial modes, we use the term of "inertial modes" to refer to 
the inertial modes with the comparable spheroidal and toroidal components of the velocity field, 
and the term of "r-modes" to refer to the inertial modes whose toroidal component of the velocity 
field is dominating the spheroidal component. 

Very recently, Lockitch &; Friedman (1998) calculated the inertial modes, assuming an 
ordering law given by 

5; -0(1), i^/~0(l), iT;, -0(1), 

5pr^O{Q?), 6p^0{n^), (5$ ~ 0(^2), a~0(J)), 

where Si and Hi are the spheroidal components and iTi/ is the toroidal component of the 
displacement vector, and Sp, Sp, and 5^ are the Eulerian perturbations of the density, 
pressure, and gravitational potential, and a is the oscillation frequency observed in an inertial 
frame. For r-modes, we may assume an ordering law given as cr ~ 0{Q,), iTi' ~ 0(1), and 

~ iJ; ~ (5/9 ~ 5p ~ (5$ ~ 0{fP). Lockitch &; Friedman (1998) showed that a number of inertial 
modes are unstable because of the gravitational radiation reaction, just as r-modes are, and 
estimated the dissipation timescales due to the gravitational radiation and some viscous processes 
for slowly rotating polytropic neutron star models, where they made use of eigenfrequencies and 
eigenfunctions for slowly rotating Maclaurin spheroids to classify the inertial modes. Notice 
that Lindblom Sz Ipser (1998) obtained analytically eigenfrequencies and eigenfunctions of the 
inertial mode for the Maclaurin spheroids (see also Bryan 1889). Lockitch & Friedman (1998) and 
Lindblom &; Ipser (1998) called the inertial modes "rotation modes" or "generalized r-modes", 
which are exactly the same as the inertial modes that will be discussed in this paper. 

In this paper, we study the inertial modes and r-modes of slowly rotating stars, by employing 
a different numerical approach to the problem (Lee & Saio 1986, Lee 1993) from the calculations 
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by, e.g., Lindblom et al. (1999), and Lockitch & Friedman (1998). Newtonian gravity, uniform 
rotation, and isentropic structure of the equilibrium are assumed for adiabatic perturbations in this 
paper. The effects of the rotational deformation of the equilibrium structure on the eigenfunctions 
and eigenfrequencies are included in the analysis, where the formulation by Lee (1993) (see also Lee 
Sz Saio 1986) has been improved so that the eigenmodes are consistent up to the order of Q^. The 
plan of this paper is as follows. In §2, we describe briefly the basic equations for the equilibrium 
configurations of slowly rotating polytropes. In §3, we present the formalism to calculate normal 
modes of slowly rotating stars. In §4, the properties of the eigenfrequencies and eigenfunctions of 
the inertial modes are discussed. In §5, we examine the stability of simple neutron star models 
against the inertial modes and the r-modes, computing their growth or damping timescales due to 
the gravitational radiation reaction and some viscous damping processes. §6 is for discussions and 
conclusions. 



2. Equilibrium State 

We consider oscillations of an uniformly rotating, isentropic, Newtonian star. The structure 
of an equilibrium state is determined by the hydrostatic equation, the Poisson equation, and the 
equation of state: 

Vip = -pV,^, (2) 

V,V* ^ = AnGp , (3) 

P = P{p) , (4) 

where ^ is the effective potential defined by 

^ = il7V2{l -P2(cos0)}, (5) 

and 0, is the angular rotation frequency of the star, constant for uniform rotation. Here we use 
spherical polar coordinates {r,6,(j)), and Pi{cos6) denotes the Legendre polynomial, and Vj is the 
covariant derivative with respect to the quantity followed. 

In this investigation, assuming slow rotation, we apply the Chandrasekhar-Milne expansion 
(see, e.g, Tassoul 1978) to equations (|2|)-(^). In this technique the effects of the centrifugal force 
and the equilibrium deformation are treated as small perturbations to a non-rotating spherically 
symmetric star. The small expansion parameter due to rotation is chosen as = i7(i?^/GM)^/^, 
where R and M are the radius and the mass of the non-rotating star, respectively. To the lowest 
order effects of the centrifugal force, the effective potential ^ can be expanded as follows: 

^(r, e) = ^oir) - 2R^n^ [Mr/R) + ^2V'2(r/i?)P2(cos 9) + 0(^2)] , (6) 

where ^o{f) is the potential of the non-rotating star. The functions ipo{x) and ip2ix) are 
determined by solving the following ordinary differential equations: 

J_d_ 

x"^ dx 



dx 



k{x)'ipo{x) + 1 , 



(7) 
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k{x)=47TGR^^, (9) 

where x = r/i? and po is the density of the non-rotating star. For simpHcity, we consider a 
sequence of slowly rotating stars whose central density is the same as that of the non-rotating star. 
From this condition, together with the regularity condition at the center of the star, the boundary 
conditions of equations (0) and (|^) at the origin are given as 

*,(0) = 0. ^M?)=o, (10) 

where i = 0, 2. From the boundary conditions at the stellar surface, the constant A2 is determined 
as ^ 

^'""6 (3^2(1) +dV'2(l)/d In x) ■ ^-^^^ 



3. Perturbation Equations 

In the perturbation analysis, we introduce a parameter a that is constant on a distorted 
effective potential surface. In practice, the parameter a is defined such that 

^{r,e) = ^o{a). (12) 

With this parameter a, the distorted equi-potential surface may be given by 

r = a{l + e{a,e,(l))}. (13) 



By using equations (g), (|l^) and ([l^), we obtain the explicit expression for the function e{a,9) up 
to the order of Cl^: 

e{a,9) = a{a) +/3{a)P2{cose), (14) 

where 

2cil7Vo(x) 

a{a) = ^ , (15) 

(3{a) = 2 • (16) 

Here, ci = {a/ Rf /{M{a)/M), and M(a) denotes the mass inside the a-constant surface. 

Hereafter, we employ the parameter a as the radial coordinate, instead of the polar radius 
coordinate r. In this coordinate system {a,6,ip), the metric tensor is written by 

ds^ = (1 + 2e)(da^ + a^dO^ + sin^ 9d(p^) + 2ae^ada^ + 2aefidade + 0(f]^) , (17) 
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where a comma denotes a partial derivative with respect to the variables followed. Note that in 
this frame the pressure, density and effective potential of a rotating star depend only on the one 
variable a, although the orthogonality of the basis vectors is lost. 

The governing equations of nonradial oscillations of a rotating star are obtained by linearizing 
the basic equations used in fluid mechanics. Since the equilibrium state is assumed to be stationary 
and axisymmetric, the perturbations may be represented by a Fourier component proportional to 
^i(at+m(fi) ^ ^]^gpg fj jg f^i^Q frequency observed in an inertial frame and m is the azimuthal quantum 
number. The continuity equation may be linearized to be 

Sp = -ViipC), (18) 

where is the Lagrangian displacement vector, and we have made use of 6v^ = i(cT + mri)^* = ito^,^ 
with CO being the oscillation frequency observed in the corotating frame of the star. The linearized 
Euler's equation is 

- LO^^i + V, f — + 6^] + A — + e^i-Vip + 2iujnC^Vnpi = , (19) 
VP/ P p 

where 99* is the rotational Killing vector, with which the 3-velocity of the equilibrium fluid of a 
rotating star is given as = Qip^. The last term on the left-hand-side of equation ( [T^ ) comes from 
the Coriolis force, and its explicit components in the metric (p!7[) are given by 



^iVjV^a = -(1 + 2e + ae,a) asin^ O^^^ + 0{n^) , (20) 

C^Vjipe = -{(1 + 2e) cos 6 + e,e sin 6 } sin + 0{n^) , (21) 

^^Vjif^ = -(1 + 2e + ae^a) a sin^ 9^" - {(1 + 2e) cos 9 + e,^ sin 9}a'^ sin 6*^^ + 0(0.^) . (22) 

(Note that Lee (1993) ignored the terms of e, e^^, and e^g in ^.-''VjLpi, as a result of which the 
eigenmodes obtained are consistent only up to the order of Q'^. By retaining these terms in the 
analysis, the eigenmodes calculated in this paper are consistent up to the order of 0^ for slow 
rotation.) The perturbed Poisson equation is given by 

ViVM$ = 4ttG6p . (23) 

For adiabatic oscillations, we have 

Sp = Tp(^^ + eA,>j (24) 



where F is the adiabatic index given by 



,(91np 

and Ai is the generalized Schwarzschild discriminant defined by 



ad 



Ai = -Vip-^Vip. (26) 
P 
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In the case of isentropic stars, the 1-form Ai vanishes exactly. 

Physically acceptable solutions of the linear differential equations are obtained by imposing 
boundary conditions at the inner and outer boundaries of the star. The inner boundary conditions 
are the regularity condition of the perturbed quantities at the stellar center. The outer boundary 
conditions at the surface of the star are Ap/p = and the continuity of the perturbed gravitational 
potential at the surface to the solution of VjV*5$ = which vanishes at infinity. 

In order to solve the system of partial differential equations given above, we employ a series 
expansion in terms of spherical harmonics to represent the angular dependence of the perturbed 
quantities. The Lagrangian displacement vector, is expanded in terms of the vector spherical 
harmonics as 

oo 

r = E aSi{a)Yrie,ip)e^''\ (27) 

l>\m\ 

e= f: U{a)Yr,e + TKa)^^Yr^e^^', (28) 

l,l'>\m\ 
1 oo 

^'' = -^ E {Hiia)Yr\^-T,ia)smeYr^,}e^'^' (29) 

l,l'>\m\ 

(Regge &: Wheeler 1957, see also Thorne 1980). The perturbed scalar quantity such as 6^ is 
expressed as 

oo 

5$= ^ 6<^i{a)Y{^{9,ip)e'''K (30) 

l>\m\ 

Substituting the perturbations into the linearized basic equations (|l^), ([l9|), (23) and (^4|), 



we obtain an infinite system of coupled linear ordinary differential equations for the expansion 
coefficients. The details are given in the Appendix. 

For numerical calculations, we truncate the infinite set of linear ordinary differential equations 
to obtain a finite set by discarding all the expansion coefficients associated with I higher than Imax- 
This truncation determines the number of the expansion coefficients kept in the spherical harmonic 
expansion of each perturbed quantity. We denote this number as kmax- Note that the number 
kmax is equivalent to the dimension of the column vectors y^, introduced in the Appendix, for 
i = 1-4. Our basic equation therefore becomes a system of 4 x /cmax-th order ordinary differential 
equations, which, together with the boundary conditions, are numerically solved as an eigenvalue 
problem of w = cr + mO, by using a Henyey type relaxation method (e.g., Unno et al . 1989, see 
also Press et al. 1992). 

For the number Aj^ax of the expansion coefficients kept in the eigenfunctions, we employ k^^^=h 
or 6, values that are found to be large enough for the inertial modes calculated in this paper. Note 
that if we consider higher radial-order modes, with many radial nodes of the eigenfunctions, than 
the modes obtained in this study, the effects of truncation can be more serious, as shown by Lee 
et al. (1992). 



- 7- 



When $7 = 0, nonradial oscillation modes are classified into two decoupled sets, called the 
"polar" parity mode (or the "spheroidal" mode) and the "axial" parity mode (or the "toroidal" 
mode). Although these modes have been traditionally called "even" and "odd" parity modes, 
respectively, we will use these terms for different meanings in this paper (see below). When $7 7^ 0, 
both parity modes are mixed, and an oscillation mode contains contributions from the polar and 
axial parity modes. 

Since the equilibrium state of a rotating star is invariant under the parity transformation 
defined by ^ — tt — ^, and is symmetric with respect to the equator, the linear perturbations have 
definite parity for that transformation. In this paper, a set of modes whose scalar perturbations 
such as are symmetric with respect to the equator is called "even" modes, while a set of modes 
whose scalar perturbations are antisymmetric with respect to the equator is called "odd" modes 
(see, e.g., Berthomiue et al. 1978). For positive integers = 1, 2, 3, • • •, we have I = \m\ + 2k — 2 
and r = Z + 1 for even modes, and / = \m\ +2k — 1 and Z' = Z — 1 for odd modes, where the 
symbols I and I' have been used to denote the spheroidal components and toroidal components 
of the displacement vector respectively (see equations (27)— (30)). Notice that in Lockitch k, 
Friedman (1998), the terms "even" and "odd" modes are used to denote the "polar-led hybrids" 
and "axial-led hybrids" modes, respectively. 



4. The Eigenvalues and Eigenfunctions of Inertial Modes 

We compute the eigenvalues and eigenfunctions of inertial modes and r-modes of rotating 
polytropic stars with accuracy up to the order of Cl^. In this study, we concentrate our attention 
to the case of isentropic stars, for which the generalized Schwarzschild discriminant Ai vanishes 
exactly and the equilibrium state is marginally stable against convection. In this case, the 
adiabatic index is given by using the polytropic index n as 

r = ^ = ^. (31) 

amp n 

Since the rotational deformation is of the order of J)^, the effects on the frequencies (oc $7) of the 
inertial modes and r-modcs appear as terms of the order of Ct^ . For sufficiently small values of f7, 
therefore, we may expand the eigenfrequency oj of the inertial modes and r-modes of a rotating 
star in powers of $7 as 

^ = ^0 + ^2^^ 0{a^) . (32) 

In this study, the expansion coefficients kq and K2 are obtained by fitting the eigenfrequencies 
computed for a mode at several small values of $7. 

For nonradial p-, /-, and 5-modes of a non-rotating or sufficiently slowly rotating star, we may 
specify the values of I and m without any ambiguity and, by counting the number of radial nodes 
of the eigenfunctions, we may order the modes without any confusion (e.g., Unno et al. 1989). But, 
to our knowledge, no well-established classification scheme for inertial modes of a rotating star 
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exists, except for the case of the Maclaurin spheroids, for which exact eigensolutions are known 
(Bryan 1889; Lindblom & Ipser 1998). As shown by Bryan (1889), eigenmodes of the Maclaurin 
spheroid can be characterized by the eigenvalue oj/Vl and the angular quantum numbers and m 
of the associated Legendre functions of the spheroidal coordinate defined on a surface form of the 
spheroid. The functional forms of the eigenfunctions in the spheroidal coordinate, which depends 
on the eigenvalue w/J^, are specified by Iq and m. Since the the spheroidal coordinate employed is 
dependent on w/f], the functional forms in the actual spatial coordinate are different for different 
eigenvalues uj/Vi for given Iq and m. In fact, Lindblom &: Ipser (1998) have shown that, for given 
Iq and m, the eigenvalue kq for inertial modes of the Maclaurin spheroid can be determined, to 
lowest order of ll, by solving the algebraic equation given by 

Jm / ^ \ ^m+1 

"^^^/o('^o/2) + (y - Ij d^^'o(^o/2) =0, (33) 
where m > is assumed. They also have shown that 

- 2 < Ko < 2 , (34) 

and that the number of the roots of equation ( [3^ ) is equal to Iq — m. 

In this paper, as a classification scheme for inertial modes of uniformly rotating polytropic 
stars with n 7^ 0, we follow the scheme used for the inertial modes of the Maclaurin spheroid, and 
we pick up the inertial modes that are similar in the mode character to those of the Maclaurin 
spheroid, labeling them with the quantum numbers Iq and m. In practice, the value of Iq is 
determined by examining the properties of the eigenfunctions, such as the number of radial 
nodes of the expansion coefficients, and the number of the dominant expansion coefficients of the 
eigenfunctions. A similar procedure was employed by Lockitch &: Friedman (1998). Notice that 
our definition of Iq is not the same as that of Lockitch & Friedman (1998), but the same as that of 
Lindblom & Ipser (1998). 

In Table 1 we tabulate the eigenvalues (kq, H2) for inertial modes and r-modes of a polytropic 
model with n = 1. As expected from the results for the Maclaurin spheroid, the value of Iq — \m\ is 
odd for odd modes and even for even modes. The modes which satisfy the condition a{a + mVt) < 
are unstable to the gravitational radiation reaction (see the next section), and are marked with an 
asterisk *. We find that a number of inertial modes are unstable to the gravitational radiation 
reaction, as first shown by Lockitch &; Friedman (1998). The modes with Iq — \m\ = 1 are r-modes 
and the frequency to lowest order in Q is given by kq = 2/{\m\ + 1), which is independent of 
the structure of the star. Only the r-modes with Iq — 1 = \m\, for which the toroidal component 
iTii associated with Y^l^^^^ is dominating, are found in isentropic stars, as suggested by Saio 
(1982). Table 1 also shows that, although r-modes are all retrograde, there are both prograde and 
retrograde inertial modes. 

In Table 2, the eigenfrequencies {ho,K2) of inertial modes and r-modes with m = 2 are 
tabulated for several values of the polytropic index n, to see the dependence of the eigenfrequencies 
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on the equation of state, where the eigenfrequencies for the case of n = are obtained by using the 
results by Lindblom Sz Ipser (1998). As shown by Table 2, the eigenvalues (K0A2) of the inertial 
modes depend on the polytropic index n, while only the eigenvalue K2 of the r-modes is dependent 
on the index n. 

Let us compare our calculations to those by Lindblom ct al (1999) and by Lockitch & 
Friedman (1998). We find that the r-modc eigenfrequencies {k,q,k,2) given in Table 1 for m =2 
and 3 are in good agreement with those given in table 1 of Lindblom et al. (1999). Note that the 
normalization applied to K2 employed by Lindblom et al. (1999) is different, by a factor 4/3, from 
that employed in this paper. We also find that the inertial mode frequencies kq given in Tables 1 
and 2 are in good agreement with those given in the table 6 of Lockitch & Friedman (1998), who 
do not give any number for K2 because of their neglect of the effects of the rotational deformation. 

In Figures 1 to 3, we show the several expansion coefficients Si{a/R), Hi{a/R) and iTi'{a/R) 
for four different m = 2 inertial modes of a polytropic model with n = 1, where the amplitude 
normalization is given at the surface of the model by iT\^m\+i — 1 even modes and iTj^j = 1 
for odd modes. In Figure 1, the first two expansion coefficients are plotted, versus the fractional 
radial coordinate defined as a/R, for the even parity inertial modes with m = 2 and Iq — |m| = 2. 
The solid curves give the expansion coefficients corresponding to the mode with = 1.100, 
while the dotted curves correspond to the mode having kq = —0.557. This figure shows that 
only the first expansion coefficients are dominating for both axial and polar components of the 
displacement vector. In Figures 2 and 3, the first three expansion coefficients of the spheroidal 
and toroidal components of the displacement vector are plotted against the fractional radial 
coordinate a/i? for the unstable even parity inertial modes with m = 2 and Iq — I'm] = 4. Figures 
2 and 3 correspond to the mode with kq = 1.520 and with kq = 0.863, respectively. The figures 
show that only the first two expansion coefficients of the eigenfunctions have dominant amplitude. 
Figures 1 to 3 may confirm that the series expansion of the eigenfunctions in terms of vector 
spherical harmonics converges quickly for the inertial modes calculated here. We note that the 
behavior of the eigenfunctions is very similar to that of the corresponding modes of the Maclaurin 
spheroid (see also Lockitch Sz Friedman 1998). If we consider the modes associated with given Iq 
and m, the number of nodes of the dominant expansion coefficients in the radial direction is the 
same for the modes with different eigenvalues kq. This means that the number of radial nodes of 
the eigenfunctions is not necessarily a good quantum number to label the modes with different 
eigenvalues kq but with the same Iq and m. 

Figures 1 to 3 also show that for inertial modes of a slowly rotating isentropic star, the toroidal 
expansion coefficients iT;/ has comparable amplitude to the spheroidal expansion coefficients Si 
and Hi in the interior of the star. It is also shown that since the eigenfunctions Si have vanishing 
amplitude at the surface, the fluid motion is almost tangential to the stellar surface. 

In Table 3, we tabulate the number of radial nodes of the dominant expansion coefficients 
of of the r-modes and inertial modes for several values of Zq — Note that the content of 
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Table 3 is the same for m = 1, 2, 3. The number of the dominant expansion coefficients of the 
eigenfunctions iTi' may be given as the maximum positive integer k that satisfies 2k — l<lQ — \m\. 
Also the number of the dominant expansion coefficients of the spheroidal components, Si and Hi, 
may be given as k for the even parity modes and as k — 1 for the odd parity modes. Although 
we cannot give any proof in general, Table 3 may show that the number of radial nodes of the 
dominant expansion coefficients iTii with the highest harmonic index is zero, and that the 
number of radial nodes of the dominant expansion coefficients iTi/ with the lowest harmonic index 
/' is given by the maximum positive integer k which satisfies 2k + 1 < Iq — \m\. 



5. Dissipation Timescales 

As mentioned in the previous section, some of the inertial modes calculated in this paper 
are unstable to the gravitational radiation reaction in the sense that the frequency a satisfies the 
condition a{a + mO) < (see, equations (|4^) and (|4^)). However, since in neutron stars there are 
some possible viscous and dissipative processes which tend to suppress the instability, to decide 
whether the rotating neutron star is really unstable due to the gravitational radiation reaction, we 
need to know that the instability is strong enough to survive the dissipative processes. 

The effects of the gravitational radiation reaction and viscous processes on the r-modes have 
already been studied by a number of authors (Lindblom et al 1998, Owen et al 1998, Andersson 
et al 1998, Kokkotas & Stergioulas 1998, Lindblom et al 1999), and they have found that for the 
r-modes the instability due to the gravitational radiation reaction strongly dominates the viscous 
damping processes considered. For inertial modes, Lockitch & Friedman (1999) suggested that 
the inertial mode instability due to the gravitational radiation survives the dissipative processes 
considered for the case of r-modes, but the instability itself is weaker than that of the r-modes. In 
this section, we reconsider these estimations using a different numerical approach to compute the 
oscillation modes of a rotating star. 

To estimate the dissipation timescales associated with viscosity and gravitational radiation 
reaction, we employ a simple method used for the analysis of the r-mode instability (Lindblom 
et al 1998, see also Ipser &: Lindblom 1991). When the dissipation timescale is sufficiently 
long compared to the oscillation period of the mode, the growth rate or the damping rate is 
approximately evaluated with the non-dissipative eigenfunctions as 



where E is the energy of the oscillation, observed in the corotating frame, given by 

(i^x , (36) 



2 



pSv'Sv* + [^ + 6<l>]6p* 



and the asterisk * denotes the complex conjugate of the indicated quantity. The time derivative of 
E is determined by the dissipation effects, for which we consider the shear viscosity, bulk viscosity 
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and gravitational radiation. 

The dissipation rate due to the shear viscosity can be calculated from 

(^)^ = -2 / vSa'^Sa*jd^x, (37) 
where (5(7 jj is the shear of the perturbation written by 

Saij = \ {Vi^Vj + Vjbvi - lgijVk6v''^ , (38) 
and the coefficient of shear viscosity for hot neutron star matter is given by 



"-^x^O'MjoT^j'f^) (39) 



(Cutler & Lindblom 1987; Sawyer 1989). The dissipation rate due to the bulk viscosity can be 
written as 

csese* d^x , (4o) 



dE 



B 



where 59 is the expansion of the perturbation defined as 



5e = vM =-#^ i^p + e^] , (41) 



Tp \ da J 

and the bulk viscosity coefficient for hot neutron star matter is 



a + mnj Vl0i^g-cm-3y \lO^K 



C = 6xl0- --— g-cm-.s- (42) 



(Cutler & Lindblom 1987; Sawyer 1989). 

The dissipation due to gravitational radiation comes from two divisible contributions, which 



are 



and 



^ = -a{a + mn)J2Nia''^mm\\ (43) 



G-D 



(7 7^ \ CO 
= -a{a + mn) ^ Nia^' \dJim\' , (44) 
"I / G-J 

where the coupling constant Ni is given by 

' c2«+i/(/-l)[(2/ + l)!!]2- ^ 

Here, the mass, SDi^^, and current, SJi^, multipole moments of the perturbation are given by 
(Thorne 1980, Lindblom et al. 1998) 



j 6pr%*"'d^x , (46) 
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and 

1 



5 Jim = \ ' / iP^^i + ^P^i) Km*d^x , (47) 



i B 

and is the magnetic type vector spherical harmonic given by 



(Thorne 1980). Notice that as indicated by equation (p7[), the gravitational radiation due to the 
current multipole moment is emitted by the "axial" parity perturbations. 



For inertial modes, since the first term in equation ( |3q ) is dominant and the energy E of 
the modes is positive definite for sufficiently small fi, the instability sets in when dE/dt > 0. As 
shown by equations (^) and (|4^), since only the energy change rate due to gravitational radiation 
can be positive and the energy change rates due to the other dissipative processes are all negative, 
the necessary condition for the instability is that there exist modes whose frequencies satisfy 



a{a + mVt) < and the right-hand-side of equations (^) and (44) becomes positive. Physically, 
this means that the rotating star becomes unstable when the gravitational radiation carries away 
negative energy. 



We may write the damping timescale of the mode to the lowest order in as follows: 
1 1 / W^kV I / T \^ firGp 



T 




+ >:-- +>:^ , (49) 



where p is the average density of the star. Here, the first, second, third and fourth terms in 
equation ( ^9[ ) are contributions from the shear viscosity, the bulk viscosity, the current multipole 
radiation and the mass multipole radiation, respectively. The expression of the timescale r given 
by equation ( ^9|) is basically the same as that given by Lockitch & Friedman (1999), who ignored 
the last term, considering that the gravitational radiation from the mass multipole moment is 
negligible for inertial modes. However, for the even parity modes, the largest components from 
the mass multipole moment and the current multipole moment may be (5D|^| and SJ^^^^i^^, 
respectively. If both of the components have the same amplitude, the radiation from the mass 
multipole moment is by a factor larger than that of the current multipole moment. Since the 
square of the mass multipole moment is by a factor Q'^ smaller than that of the current multipole 
moment because of the ordering law 6p ~ 0(^1^) and 6v^ oc o;^* ~ 0(fi) given by equation (1), the 
radiation reaction from the mass multipole is the same order as that from the current multipole 
for the even parity inertial modes. For consistency, therefore we must include the mass multipole 
moment as well as the current multipole moment to estimate the energy change rate. 



In Table 4, we tabulate the timescales in the unit of second for various dissipative processes 
for a polytropic model with the index n = 1, where the radius and the mass at = are chosen to 
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he R = 12.57km and Af = 1.4Mq, respectively. As shown by Table 4, for inertial modes the mass 
multipole radiation dominates the current multipole radiation for both the even and odd parity 
modes. Note that for r-modes with Iq — \m\ = 1, only the current multipole moment is relevant 
to the gravitational radiation. We also find that for even parity inertial modes, the fastest growth 
time due to the mass multipole radiation is of the order of 10^ seconds for the parameters 10^ K 
and Q = ^'KGp. 

Most of the timescales calculated in this paper are in good agreement with those obtained 
by Lockitch & Friedman (1998) and Lindblom et al. (1999), and the relative difference between 
the calculations are at most several percents. However, there are a few discrepancies between 
our calculation and the calculations by Lockitch & Friedman (1998) and Lindblom et al. (1999). 
We note that the most dominating current multipole 5J/m is that associated with Imax and the 
moment bJim with I < Imax is less important than SJi^^^m, where Imax is the largest value of I 
associated with the dominating expansion coefficients of 6v°'. The result that for the odd parity 
inertial modes the multipole moment 6 Jim with / < Imax does not vanish may contradict the result 
obtained by Lockitch & Friedman (1998), who suggested that for the odd parity inertial modes all 
the current multipole moments 5 Jim with / < Imax vanish (or nearly vanish). The origin of these 
discrepancies is not clear at the moment. 



As seen from equation (|49|), the total dissipation timescale r is given as a function of the 
angular velocity and the temperature T of a star. For a given temperature T, we can define the 
critical angular velocity flc such that 

for < Jlc < (vrGp)^/2 _ (50) 



r(r!c,T) 

With this critical f^c, we may say that a rotating neutron star is unstable due to the gravitational 
radiation when > because the gravitational radiation reaction dominates the viscous 
damping processes. In Figure |^, the critical angular rotation velocities ilc ior two inertial modes 
and an r-mode are plotted against the temperature of the star. From this figure we can see that 
the instability of the m = 2 even parity inertial mode with kq = 1.100 is strong, although this 
instability is weaker than that of the r-mode. At about T = lOpK the critical angular velocity of 
this mode is about 25% of (vrGp)-*^/^. 



6. Discussion and Conclusions 

In this paper, we have investigated the properties of inertial mode and r-mode oscillations in 
rotating isentropic Newtonian stars, using the slow rotation approximation. By taking account 
of the effects of the rotational deformation of the equilibrium state, we have evaluated the 
eigenfrequencies of the inertial modes and the r-modes with an accuracy up to the order of 0^. 
We have also estimated the dissipation timescales due to the gravitational radiation and the 
viscosity for a simple neutron star model. We show that the inertial modes emit gravitational 



-14- 



radiation mainly by the mass multipole rather than the current multipole. It is also found that 
the instability due to the gravitational radiation reaction is strong for the most unstable incrtial 
mode, although the instability associated with the inertial modes is not as strong as that with the 
r-modes. 

In spite of the recent improvements in our understanding about the instability of the inertial 
modes and the r-modes, it seems that the fundamental properties of these modes have not yet 
been sufficiently understood, considering that most of the previous investigations of the inertial 
modes and the r-modes, including the present paper, are restricted to the case of uniformly and 
slowly rotating, isentropic, Newtonian stars. In this sense, we do not have any clear understanding 
about the inertial modes and the r-modes of, for example, differentially and rapidly rotating, 
non-isentropic, relativistic stars. As a first step to relax these restrictions, several studies have 
already appeared. For example, Andersson, Lockitch &; Friedman (1999) suggest a possibility that 
the assumption of the purely axial mode oscillation leads to inconsistent radial behavior for the 
case of isentropic relativistic stars. Kojima (1998) and Kojima Sz Hosonuma (1999) also argue 
that the r-mode spectrum is continuous for relativistic stars. Therefore it is highly desirable to 
investigate the properties of the inertial modes and the r-modes in the more general case than that 
investigated previously and in this paper, to conclude that the instability does have an importance 
in astrophysics. 

We would like to thank Prof. H. Saio for useful comments. S.Y. would like to thank Prof. Y. 
Eriguchi, Prof. T. Futamase and Dr. S. Yoshida for discussions. S.Y. was supported by Research 
Fellowship of the Japan Society for the Promotion of Science for Young Scientists. 



The derivation of the governing equations of nonradial oscillations of rotating stars is almost 
the same as that given by Lee & Saio (1986) and Lee (1993), except that the governing equations 
given in this Appendix are formulated to calculate the eigenfrequencies and eigenfunctions correct 
up to the order of $7^. 

We introduce column vectors yi, y2, ys, y4, h, and t, whose components are defined by 



A. Basic equations for nonradial oscillations of slow^ly rotating stars 



yi,k = Si{a) 



(Al) 




(A2) 



y3,k 




(A3) 



(A4) 
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h,k = Hi{a) , (A5) 

and 

t,k = Tv{a), (A6) 

where 

,= -'-fia). (A7) 

p da 

Here I = \m\ + 2k — 2 and I' = I + 1 for "even" modes, and / = \m\ + 2k — 1 and I' = I — 1 for 
"odd" modes, and in both the cases k denote the indices of the vectors, and = 1, 2, 3, . . .. 

In vector notation, equations of the adiabatic nonradial pulsation for a sfowly rotating star 
are written as foHows: 

The perturbed mass conservation law ( |l8|) reduces to 

+ 3??(/?) ^0 h + 3mi?(/3) Qo «t . (A8) 
The a component of the perturbed Euler's equation (|l9| ) reduces to 

a'^^ - iciu? + aAa)yi- [l- aAa-U)y2- aAayz + 2mciu}Oh. + 2ciu:ClCQit 
da 

= ci(:D2[2{7?(a)l +r/(^)A}yi -3/?eoh-3m^Qo«t 
-mv{{a + r?(a))l + {(3 + h 
-v{{a + 7?(a))Co + {(3 + ??(/?) )^:o} ^t] . (A9) 

The perturbed Poisson equation (p^) reduces to 

a^-(i-f/)y3-y4 = o, (Aio) 



and 



Here, 



dyA V f V \ 

+ aAaUyi-U-y2-{^o-U-l\Ti + Uyi 

= 2{ry(a)l + r?(/3)A} a ^ + {F(a)l + A - 6/3(A + ^o)} y4 

da 

- 2{al + /3A}Aoy3. (All) 



C/ = ^^, ^ = -^, (A12) 

a In a a In a 

da da 

via) = 3a + o— , 7?(a) = a + a— , (A13) 

da da 

F(a)=2U'qia)-a—+a—[a—], (A14) 

aa da \ da J 
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uj = uj / [GM / R^Y^"^ is the frequency in the unit of the Kepler frequency, and v = 20./ uj . 
The 9 and (p components of the perturbed Euler's equation ( |l9| ) reduce to 

Loh-Miit = ^y2 + {6 + MiLr'K}yi, (A15) 

Li it - Mo h = -K yi , (A16) 



where 



Lo = (1 + 2a)Lo + /3Ao ^ Wo - 2mi//3Ao ^ (6 A + 1) , (A17) 

Li = (1 + 2a)Li + /JA^-^Wi - 2mv(3A.l^ (6 ^ + 1) , (A18) 

Mq = (1 + 2a - 2/3)Mo - 6m/3Af ^Qi + Au^A.^^ (PiAq + 3£:i + Ci) , (A19) 

Ml = (1 + 2a - 2/3)Mi - Qm[3K^^Qo + 4z^/3Ao ^ (PoAi + S^q + Co) , (A20) 

K = (1 + a + 7?(a)) z^K - Sm/^A^^Qi + + r/(/3)) A^^ (4Pi + ,?i - 2Qi) , (A21) 



and 



6 = A^i 



mv[{l + a + r/(a)) ! + (/? + r]{f3)) Aq] - 3/3 (2^ + Bq) 



MiLj^^K; (A22) 



Wo = 2AAo + 6i3o , Hi = 2^1 Ai + 6i3i . (A23) 

The quantities Qo, Qi, Co, Ci, K, Lo, Li, Ao, Ai, Mo, Mi, Aq, Ai, Bq, Bi, Vq, Vi, Sq, and 
£i are matrices written as fohows: 



For even modes, 



(Qo)i,i = Jl+l , (Qo)i+l,i = Jl+2 ' 

(Qi)j,i = Ji^i , (Qi)j,i+i = Ji+2 ' 

(Co)i,i = -{I + 2) J^+i , (Co)i+i,i = (/ + l)J/+2 ) 

(Ci)i,i = IJl^li , (Ci)i,i+i = + 3)J/+2 ; 

(K),, = ^, (K),,+i = -^, 



/(/ + !)' ' ''''' (/ + !)(/ + 2)' 

(Ao)i,i = l{l + 1) , (Ai)i,i = {l + + 2) , 

(Mo),. = Jr+, , (Mo).,+i = z.[±| J,™2 , 

(Ml),, = Jll, , (Mi),+i, = Jj±^ J,-2 , 

1 3 2 2 3 

(A)i,i = ~2 ~^ 2 (-^i+i) } ' (A)i+i,i = (A)i,j+i = 2 ) 
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(^o)v = -{1 + i){jr? , (SoUi, = ijr+iJr+2 , 
{Bo)i,i+i = -{I + 3) j,-2 , (^i)v = {1+ mji^2? -{1+ 'i){jr+i? , 

m\2 I /' rm \2 i ( jm \2i 



jm 



(^o)i+l,i = 4+2 
= Jill 
iPl)i,i+l = Jl+2 

{£o)i+l,i = Jl+2 



l + l{ijrr+iJ'iiir+iJn2y} 



jm jm jm 

0)i,i+l — 2 '^l+l'^l+2'Jl+3 ' 



l + l{{Jlllf + iJn2f + iJn3?} 



, in 

, (^o)i+2 



rm Tm rm 
2 •^«+2-'i+3'^i+4 



/-T-) \ rm Tm jm 

/T-) \ rm jm jm 

, \'^\)i,i+2 — 2 <^i+2<^Z+3<^Z+4 ) 



i (/ + 2) + ^ {(/ + i)(j,72)' -{1+ 2){{jrf + (J;'?!)')} 



4 + 1) + ^ {(^ + l)(('^/+2)' + (4^3)') -(1 + 2)(J^Tl)'} 



3 3 

(^o)i,i+l = "2 + '^)Jl+lJl+2Jl+3 ' (^o)i+2,i = 2 l)4+2'^Z+34+4 > 



(^l) 



1 



^ + ^ + ijn2f) -{1+ wr?} 



i,i+l 



= J 



1+2 



2 a + 3) + 2 {(^ + 2)(J^?3)' - a + 3)((J,7i)' + iJr+2Y)} 



3 3 

= 2 ^Jl+lJl+2Jl+3 ) (^l)«,i+2 = — 2 ^)Jl+2Jl+3Jl+4 ' 



where Z = |m| + 2i — 2, z = 1, 2, 3, . . ., and 



{I + m){l — m) 



{21 - 1){21 + 1) 



(A25) 



For odd modes, 



Tm 
•^l+l 



jm 

"^l+l ' 



{Qo)i,i = Jr 5 (Qo)i,i+l 

= — (Z + 1) J/" , (Ci)j_|_i^j 
(K)i+i,, 



7 T-m 



jm 

(K)i,i = - — , 



jm 

■''l+l 

Z + 1 ' 



(Ao)i,i = l{l + 1) 
l + l 



{A,)i,i = lil-l), 
I 



1^22 



1 



= {i- i)jrjiii > = -{i + 2)^r^Hi ' 

m \2 I / Tm\2 , / 7m \2 



(^o)i,i+l = Ji 



m 



{V 



— 7™ 
l)i+l,i — "Jl+l 

{£0) 



-l + l{{Ji^ir+{jrr+{jnir} 

-l + l{{Ji^i? + {jrf + iJr+if} 
-l + ^iijrf + iJiiif + iJwf} 



(2^0) 



i+l,i ~ 2 ' '^l+i'^l+2 ' 



fTt \ I'm m 

\'-'Q}i,i+2 — 2 'Jl+l'Jl+2'Jl+Z 

1. 

/-n ^ 7™ 7™ 7™ 



i+2,i 



jm 

i — 'Jl 



^(^ - 1) + ^ {(^ - i)((^r)' + (Jiiif) - KJi^iT} 



2 

rm \2 



jm jm Jin 
<-';_i_i 'Jij-O'Ji 



l+l'^l+2'^l+3 



{£0)1,1+1 — Jl+1 



1(1 + 2) + ^ {-(/ + 2){{jn' + {jr+if) + i){jr+2?} 



3 3 

{£o)i+l,i ^ ~ ^'^'^l^'^l+^'^W ' {£0)1,1+2 = "2 + 4)4+l<^/+24?3 ' 



(^1) 



1,1 — 'Jl 



1(1 + 1) + I {/( J,7i)2 - (Z + l)((J,'!!i)2 + (JD^)} 



i+l,i 



l+l 



{K{jr+i? + {jr+2)^) wr?} 



3 3 

{£l)i,i+l = "2 + ^)JrJl+lJl+2 I (^l)i+2,j = 2 ^^1+1^1+2^1+3 ' 
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Eliminating h and it from equations ( [A8D and ( [ADD by using equations ( A15| ) and ( |A16 ), 
equations ([A8|)-( |ATi] ) reduce to a set of first-order linear ordinary differential equations for yi, y2, 
y3 and as follows: 



da 



V 



3 1 + yi + 



^12 V, 



V 

1 ) y2 + ^ ya 



(A27) 



dy2_ 

da 



{(ciLj2 + aAa) 1 + ciLD^^^ai} yi + {{I - aAa - U)l + ^22} y2 + aAa ys , (A28) 

(A29) 



a-r- = (1 - L/jya + y4 
da 



where 



a^ = -aAaUJ-^yi + C/ ^ J- V2 + J"' |ao - C/ ^1 - 2(al + /3A)Ao} ya 

+ J-^{-Ul + F(a)l + F{P)Ao - 6/?(A + ^o)} , 

= {1 + 3i?(/3)p} WO - a - a - 3mi?(/3)QoLr'K , 

^12 = {l + 3i?(/3)p}w, 

+ i/{(l + a + 7?(a))Co + (/3 + r/(/3))^:o} + 3m/3Qo L^^K, 



7^ 



W = Ao(Lo-MiLriMo)"S 
J = l-2{?7(a)l + 7?(/3)A}■ 
The surface boundary conditions are 

yi - y2 + ys = , 



and 



(A30) 

(A31) 
(A32) 

(A33) 

(A34) 
(A35) 



j^22 = -nw, 

p= (i3o + mQoLr'Mo)Ao\ 
mi/{(l + a + r/(a))l + (/? + ??(/?)Mo} + 3/3^0 
z.{(l + a + ??(a))Co + (/3 + 7?(/3))£:o + 3m/3Qo Lr'^oAo ' , (A36) 



(A37) 
(A38) 

(A39) 
(A40) 



A2{(1 + a)l + PAo} + {al + /3A}Ao ys + {A2(al + pAo) + l}y4 = , 

where A2 is a matrix whose elements are given by 

{A2)u = (/ + !). (A41) 

Here / = \m\ + 2i — 2 for "even" modes, / = \m\ + 2i — 1 for "odd" modes. The inner boundary 
conditions at the stellar center are the regularity conditions of the eigenfunctions. 
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Fig. 1. — The first two expansion coefficients Si{a) [panel a], Hi{a) [panel b], iTi{a) [panel c] for 
the even parity inertial modes with m = 2 and lQ — \m\ = 2 arc plotted against a/R for a polytropic 
model with the index n = 1, where we assume = 0.01. The solid curves give the expansion 
coefficients corresponding to the mode with kq = 1.100, while the dotted curves correspond to 
the mode having kq = —0.557. The eigenfunctions are normalized so that iTz{a = R) = 1. 
Attached labels Z = 2, 3, 4, 5, • • • denote the spherical harmonic indices associated with the expansion 
coefficients. 

Fig. 2. — The first three expansion coefficients Si{a) [panel a], Hi{a) [panel b], iTi{a) [panel c] for 
the m = 2 even parity inertial mode with kq = 1.520 and Iq — |m| = 4 are plotted against a/R for a 
polytropic model with the index n = 1, where Q. = 0.01. The eigenfunctions are normalized so that 
iT^{a = R) = 1. Attached labels Z = 2, 3, 4, 5, • • • denote the spherical harmonic indices associated 
with the expansion coefficients. 

Fig. 3. — The same as Figure 2, but for the mode having eigenvalue kq = 0.863. 

Fig. 4. — Critical angular velocities flc for the r-mode and two inertial modes are plotted against 
the temperature of a neutron star. The curve labeled "r-mode" is for the m = 2 r-mode with 
Ko = 0.667. The curve labeled "even" is for the m = 2 even parity inertial mode with kq = 1.100. 
The curve labeled "odd" is for the m = 2 odd parity inertial mode with = 0.517. Here, the 
critical angular velocity is normalized by using (irGp) and the unit of the temperature is Kelvin. 
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Table 1. Eigenvalues {kq, K2)^oI inertial modes and r-modes for the p = Kp^ Polytrope. 



lo - \m\^ 


m = 


= 1 


m = 


= 2 


m = 


= 3 


1" 


( 1.00000, ■ 


-0.00018) 


( 0.66667, 


0.39827)* 


( 0.50000, 


0.42727)* 


2 


(-0.40137, 


0.20020) 


(-0.55659, 


0.12248) 


(-0.63164, 


0.05993) 




( 1.41300, 


0.23767) 


( 1.10003, 


0.46560)* 


( 0.90491, 


0.54494)* 


3 


(-1.03238, ■ 


-0.21090) 


(-1.02588, - 


-0.22877) 


(-1.01487, ■ 


-0.23750) 




( 0.69059, 


0.36221)* 


( 0.51734, 


0.37950)* 


( 0.41265, 


0.35766)* 




( 1.61373, 


0.26999) 


( 1.35778, 


0.45492)* 


( 1.17674, 


0.55092)* 


4 


(-1.31227, ■ 


-0.32577) 


(-1.27289, - 


-0.33970) 


(-1.23863, ■ 


-0.34435) 




(-0.17879, 


0.01402) 


(-0.27533, 


0.00051) 


(-0.33327, • 


-0.01537) 




( 1.05153, 


0.45459) 


( 0.86295, 


0.50140)* 


( 0.73430, 


0.50562)* 




( 1.72626, 


0.24391) 


( 1.51957, 


0.41210)* 


( 1.36056, 


0.51593)* 



^Q{ko + K2^'^) = w is the mode frequency in the corotating frame up to the order 
of The modes marked with an asterisk * satisfy the condition a{a + mQ) < 
to the lowest order in O, and are considered to be unstable to the gravitational 
radiation reaction in the absence of viscous dissipation. 

''Our definition of the angular quantum number lo is not the same as that of 
Lockitch & Friedman (1998), but the same as that of Lindblom & Ipser (1998). 
Eigenmodes having odd (even) values of lo — \m\ are odd (even) parity modes. 

"These modes are r-modes. To the lowest order in fi, kq = 2/(m-|-l) (Papalouizou 
& Prigle 1978) and is independent of the equation of state. 
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Table 2. Eigenvalues (kq, K2Y of m = 2 inertial modes and r-modes for several values of the 

polytropic index n. 



It) - ' 


;) = 


= 


;) = 




;) = 




/; = 


: 2 


1" 


( 0.66667, 


0.76543)* 


( 0.66667, 


0.39827)* 


( 0.66667, 


0.28248)* 


( 0.66667, 


0.19701)* 


2 


(-0.23193, 


0.02930) 


(-0.55659, 


0.12248) 


(-0.69646, 


0.13213) 


(-0.82771, 


0.12773) 




( 1.23193, 


0.88145)* 


( 1.10003, 


0.46560)* 


( 1.06259, 


0.33590)* 


( 1.03438, 


0.24322)* 


3 


(-0.76334,- 


-0.55320) 


(-1.02588,- 


0.22877) 


(-1.12777,- 


0.15318) 


(-1.21712,- 


0.11767) 




( 0.46690, 


0.62847)* 


( 0.51734, 


0.37950)* 


( 0.53563, 


0.29081)* 


( 0.55163, 


0.22228)* 




( 1.49644, 


0.77599)* 


( 1.35778, 


0.45492)* 


( 1.31003, 


0.34136)* 


( 1.27054, 


0.25749)* 


4 


(-1.09257,- 


-0.77973) 


(-1.27289,- 


0.33970) 


(-1.34194,- 


0.24015) 


(-1.40134,- 


0.18109) 




(-0.10179, 


0.00439) 


(-0.27533, 


0.00051) 


(-0.36424,- 


0.01189) 


(-0.45689,- 


0.02235) 




( 0.88425, 


0.89693)* 


( 0.86295, 


0.50140)* 


( 0.85863, 


0.38371)* 


( 0.85688, 


0.29910)* 




( 1.64344, 


0.63981)* 


( 1.51957, 


0.41210)* 


( 1.47219, 


0.32232)* 


( 1.43071, 


0.25292)* 



^fi(K() + K2n^) = a; is the mode frequency in the corotating frame up to the order of O"*. The modes 
marked with an asterisk * satisfy the condition a{a+mQ,) < to the lowest order in and are considered 
to be unstable to the gravitational radiation reaction the absence of viscous dissipation. 

''Our definition of the angular quantum number lo is not the same as that of Lockitch & Friedman 
(1998), but the same as that of Lindblom & Ipser (1998). Eigenmodes having odd (even) values of Zq — 
are odd (even) parity modes. 

°These modes are r-modes. To the lowest order in f2, kq = 2/{m + 1) (Papalouizou & Prigle 1978) 
and is independent of the equation of state. 
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Table 3. Node number of the dominating expansion coefficients of 6v'^ for r-modes and inertial 

modes. 



lo - \m\'' 1 


2 


3 


4 


node number oi S\m\ ^ • • ' 


1 




2 


node number of 




1 




node number of S\rn\+2 






1 


node number of Si^i+s 








node number of 


1 




2 


node number of i^iml+i 




1 




node number of -H"|^|+2 






1 


node number of i/|^|+3 








node number of iT|„| 




1 




node number of iT^rn\+i 







1 


node number of iT\m\+2 









node number of ir|^|_|_3 










'^Our definition of angular quantum number lo 
is not the same as that of Lockitch and Friedman 
(1998), but that of Lindblom and Ipser (1998). 
Eigenmodes having odd (even) values of lo — \m\ 
are corresponding to odd (even) parity modes. 

''For expansion coefficients Si, we include the 
node at the stellar surface in the count of nodes. 
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Table 4. Dissipative timescales of inertial modes'^ and r-modes for n 

T = 10^ K and Q = 



1 polytropic model at 



m 


lo - m 




rs(s) 


Ts(s) 


f2(s)'^ 


Ms) 


f4(s) 


^5(s) 


1 


3 


0.691 


5.89 X 10^ 


9.23 X 10^ 


-2.46 X lO^t 


-1.27 X 10® 






2 




0.667 


2.03 X 10" 


2.50 X 10® 


-3.31 X 10° 


-3.49 X lO^t 








2 


1.100 


3.35 X 10^ 


1.23 X 10® 


-1.71 X lO^t 


-3.43 X 10" 








3 


0.517 


6.47 X 10^ 


6.18 X 10^ 


-1.31 X 10® 


-8.39 X lO^'t 


-1.88 X 10^ 








1.358 


4.13 X 10^ 


7.14 X 10^ 


-7.12 X 10^ 


-1.69 X lO^t 


-1.63 X 10** 






4 


0.863 


1.94 X 10^ 


4.92 X 10^ 


-1.63 X lO^t 


-1.12 X 10^ 










1.520 


4.82 X 10^ 


4.87 X 10^ 


-3.35 X lO^t 


-1.90 X 10" 






3 


1" 


0.500 


6.63 X 10^" 


1.43 X 10® 




-3.17 X 10^ 


-1.88 X lO^t 






2 


0.905 


1.88 X 10^ 


9.41 X 10^ 




-8.62 X lO^t 


-2.77 X 10" 






3 


0.413 


6.97 X 10^ 


4.78 X 10^ 




-1.86 X 10^" 


-5.30 X lO-'^t 


-4.07 X 10® 






1.177 


2.49 X 10^ 


6.11 X 10^ 




-1.65 X 10" 


-5.11 X lO^t 


-6.11 X 10'^ 




4 


0.734 


1.29 X 10^ 


4.10 X lO'^ 




-2.00 X 10®t 


-1.16 X 10' 








1.361 


3.06 X 10^ 


4.18 X lO'^ 




-3.45 X lO'^t 


-1.75 X 10^° 





''We present dissipative timescales only for tliose modes that are unstable to gravitational radiation reaction. 

''We present dissipative timescales fi due to the gravitational radiation reaction only for those of the dominant 
multipole moments, where for modes with odd Iq T2k+i ~ TD.2k+i and T2k ~ '?j,2fe, and for modes with even lo 
■f2fe+i = ^j,2fe+i and f2k ~ TD,2k- Here fe = 1, 2, 3, . . .. Dissipative timescales due to the mass multipole radiation are 
marked with a dagger f . 

'^This mode is the r-mode. Since the expansion of perturbations vanishes to the lowest order in 57 for r-modes, the 
definition of the bulk viscosity timescale tb for this mode is not the same as that of equation (^9|) , but that of Lindblom 
et al. (1999). 



